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Abstract 

We study the limit behaviour of solutions of 9tit — Ait + fe(|a;|) \u\^~ u = inR^x(0,T) 
with initial data kSa when fc ^ oo, where ft is a positive nondecreasing function and p > 1. 
If h{r) = r*^, (/3 > — 2) we prove that the limit function Uoo is an explicit very singular 
solution. If liminfr^o J^'^ ln(l//i(r-)) > 0, Uoc has a persistent singularity at (0, t) {t > 0). 
If r ln{l / h{r)) dr < oo, Uoo has a pointwise singularity localized at (0,0). 
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1 Introduction 

Consider 

dtu- An + h{x)\u\''~^u = in Qt := x (0, T), (1.1) 
with p > 1 and h is a nonnegative measurable function defined in M.^ . It is well known that if 

h{x) E'P {x, t)dx dt <oo, (1.2) 

r 

where E{x,t) — (47rt)~^/^e~l^l^/'** is the heat kernel, then, for any fc > there exists a unique 
solution u — Uk to (1.1 ) satisfying initial condition 

u{.,0) = k6o (1.3) 

in the sense of measures in R^. Furthermore the mapping fc i-^ is increasing. If it assumed that 
h is positive essentially locally bounded from above and from below in \ {0}, then the set {uk} 
is also bounded in the Clg^{Qrp \ {0 x (0, oo)})-topology. Thus there exist Uqo ■— hm^^oo ""fe and 
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Uoo is a solution of (1.1 ) in Qt \ {0 x (0, oo)}. Furthermore Uoo is continuous in Qj, \ {0 x [0, oo)} 
and vanishes on \ {0} x {0}. We shall prove that only two situations can occur: 

(i) Either Uoo{0,t) is finite for every t > and Uoo is a solution of (1.1 ) in Qt- Such a solution 

which has a pointwisc singularity at (0, 0) is called a very singular solution (abr. V.S.S.) 

(ii) Or Woo(0, t) = oo for every t > and Mqo is a solution of (1.1 ) in Qt \ {0 x (0, cxd)} only. Such 
a solution with a persistent singularity is called a razor blade (abr. R. B.). 

In the well-known article [4], Brezis, Peletier and Terman proved in 1985 that Uoo is a V.S.S. , if 

h{x) = 1. Furthermore they showed that Uoc{x,t) — t"^/^P~^'> f{x/ \/t) for {x,t) S Qt where / is 
the unique positive (and radial) solution of the problem 

1 -A/-i,.V/-^/ + |/rV = inM- ^^^^ 

I limi^Hoc \'n?'^-^^f{n) = 0. 

Their proof of existence and uniqueness relied on shooting method in ordinary differential equations 
(abr. O.D.E.). The already mentioned self-similar very singular solutions of the problem (1.4 ) 
was discovered independently in [6] too. Later on, a new proof of existence, has been given by 
Escobedo and Kavian [8] by a variational method in a weighted Sobolev space. More precisely 
they proved that the following functional 

V - J{v) = \l (\Vv\^ - -^^v^ + -^W^"] K{v)dv (1-5) 

achieves a nontrivial minimum in _ff^(M^), where K(ri) = e'''' 

In this article we first study equation (1.1 ) when h{x) = \x\^ {(3 € M). Looking for self-similar 
solutions under the form u{x,t) = t~^'^~^^^^'^^P~^^ f{x/\/i), we are led to 

-Af-^rj.Vf-^^^f+\vf\fr'f = inR^ 

f e HUM^) n Lf+^(M^; \vfdv) n CHR^ \ {0}) (1-6) 

lim|^l^^|,7|(2+/3)/(f-i)/(,,) = 0, 
and the associated functional 

^ ^ ^^^^ = I L ('^^'' - + (1-^) 

We prove the following 

Theorem A /- Assume f3 < N{p — 1) — 2; then there exists no nonzero solution to (1.6 ). 
II- Assume /3 > N{p — 1) — 2/ then there exists a unique positive solution f* to (1.6 ). 

One of the key arguments in the study of isolated singularities of (1.1 ) is the following a priori 

estimate 

- (^+|-^|2)(2+/3)/2(p-l) V(x,t)eQT (1.8) 

valid for any p> 1 and (3 > —2. The remarkable aspect of this proof is that it is based upon the 
auxiliary construction of the maximal solution of (1.1 ) under a selfsimilar form. Next we give two 
proofs of II, one based upon scaling transformations and asymptotic analysis of O.D.E., combining 
ideas from [4], [5] and [10], and the second based on variational methods, extending some ideas 
from [8] and valid in a more general context. As a consequence we prove 
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Theorem B Assume P > N(p - 1) - 2, then u^ (x, f) = f-'-^+^^^/^'-P^^'^ f*{x/Vt). 

It must be noticed that, if /3 < N{p — 1) — 2, Uk does not exist, and more precisely, the isolated 
singularities of solutions of (1.1 ) are removable. 

Next we consider the case of more degenerate potentials h{x): 

^ ^0 as |x| ^ Va > 0. (1.9) 

In the set of sTic;h potentials we find the borderline which separates the above mentioned two 
possibilities (i) — (V.S.S.) and (ii) — (R.B). Remark that in the case of flat potentials like (1.9), 
the corresponding solution Uoo{x,t) does not have self-similar structure and we haveto find some 
alternative techniques for the study of the structure of Uoo- Main results of the paper are the 

following two statements. 

Theorem C (sufficient condition of V.S.S. solution) Assume that the function h is continuous 
and positive in M.^ \ {0} and verifies the following flatness condition 



\x\^ In [j^j < ^ h{x) > e-^(l^l)/l^l Va; e M^, (1.10) 

where the function uj >0 is nondecreasing, satisfies the following Dini-like condition 

u)(s)ds /, 

/ ^^<oo, 1.11 
Jo s 

and the additional technical condition 

su' {s) < {2 - ao)u!{s) nearO, (1.12) 

for some ao € (0,2). Then Uoo{x,t) < oo for any {x,t) G Qt- Furthermore there exists positive 
constants Ci (i = 1,2,'6), depending only on N, ao and p, such that 

[ u^(a;,t)da; < Ciiexp [C2 ($"^ (Cai))"^] Vt > 0, (1.13) 

where is the inverse function of 

$(r) := r 
Jo 



Notice that (1.11 )-(1.12 ) is satisfied if h{x) > Ce"!^!" ^ for some 6» > 0. 

Theorem D (sufficient condition of R.B. solution) Assume h is continuous and positive in 
\ {0} and satisfies 

liminf |a;pln ( ) > <^ 3uo = const > : h{x) < exp f ) . (1.14) 

x^o \h{x)J V fIV 

Then Uoo(0,f) = 00 for any t > 0, and t Uca{x,t) is increasing. If we denote U{x) = 
limt^oo 'Woo(a^)i)) then U is the minimal large solution of 

-Au + h{x)uP = mM^\{0}, (1.15) 
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i.e. the smallest solution of (1.15 ) which satisfies 

/ u{x)dx = 00 Ve > 0. (1.16) 

Theorem C is proved by some new version of local energy method. A similar variant of this 
method was used in [1] for the study of extinction properties of solutions of nonstationary diffusion- 
absorption equations. 

Theorem D is obtained by constructing local appropriate sub-solutions. The monotonicity and 
the limit property of characteristic of razor blades solutions [16]. 

A natural question which remains unsolved is to characterize Uoo if the potential h{x) satisfies 



h(x) K, exp ( — - 
where a;(s) — > as s — > and 

uj(s)ds 
= 00. 

s 

This article is the natural continuation of [12], [14] where (1.1 ) is replaced by 

dtu- \u + h{t)\uf~'^u = Q inQr. (1.17) 

In equation (1.17 ), the function h G C([0,T]) is positive in (0,T] and vanishes only at t = 0. In 
the particular case h{t) = (/3 > 0), Uh exists if and only if 1< p < 1 + 2(1 + /3) /N, and ?i2c IS an 
explicit very singular solution. If h{t) > e^'^^*^/* where ui is positive, nondecreasing and satisfies 



^ \/u(s)ds 

= 00, 



/o s 

then Moo has a pointwise singularity at (0, 0). If the degeneracy of h is stronger, namely 

lim mi tin h(t) > —00, 

it is proved that the singularity of Uk propagates along the axis i = 0; at end, Uoo is nothing else 
than the (explicit) maximal solution ^{t) of the O.D.E. 

^' + h{t)W = in (0,00). (1.18) 

A very general and probably difficult open problem generalizing (1.1 ) and (1.17 ) is to study 
the propagation phenomenon of singularities starting from (0,0) when (1.1 ) is replaced by 

dtu- Au + h{x,t)\uf~'^ u = in Qt, (1.19) 

where h € C{Qrp) is nonnegative and vanishes only on a curve V C Qx starting from (0,0). It is 
expected that two types of phenomena should occur: 

(i) either u^o has a pointwise singularity at (0,0), 

(ii) or Uoo is singular along F or a connected part of F containing (0, 0). 

It is natural to conjecture that the order of degeneracy should be measured in terms of the 
parabolic distance to F and of the slope of F in the space x R. This could serve as a starting 
model for nonlinear heat propagation in inhomogeneous fissured media. 

Our paper is organized as follows: 1 Introduction - 2 The power case - 3 Pointwise singularities - 

4 Existence of razor blades. 

Acknowledgements The authors have been supported by INTAS grant Ref. No : 05-1000008- 
7921. 
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2 The power case 

In this section we assume that h{x) = {x]'^ with /? e M and the equation under consideration is the 
following 

dtu - Au + \xf \uf~'^ u = inQr :=K^ X (0,T) (2.1) 

with p > 1. By a solution we mean a function u € C'^''^{Qt)- Let E{x,t) = (47rf)-^/2e-l^l^/^* be 
the heat kernel in Qt and E[(^] the heat potential of a function (or measure) (j) defined by 

E[^](-'*)=(i;^X/-'^-^'^^"^(2/)'^^- (2-2) 

If there holds 

EP{x,t)\xfdxdt<oo, (2.3) 



Qt 



it is easy to prove (see [12, Prop 1.2], and [18, Th 6.12]), that for any fc e K, there exists a unique 
function u = Uk e L^{Br x (0,T)) n Lp{Br x (0,T); \xfdx) such that 

jj (^-udtC-uAC+\x\'^\uf~^uC^dxdt = kC{0,0), (2.4) 

for any ( G Cq'^(]R^ x [0,T)). By the maximum principle fc i— > Ufe is increasing. Next, it is 
straightforward that (2.3 ) is fulfilled as soon as 

f3 > max{N{p - 1) - 2; -N} (2.5) 
2.1 The a priori estimate and the mctximal solution 

In order to prove an a priori estimate, we introduce the auxiliary A'' dimensional equation in the 
variable r] = x/Vi 

- A/ - ^rj.Vf - 7/ + l/l^"' / = 0, (2.6) 

where7= (2+/3)/2(p- 1). 

Proposition 2.1 Let a > and /? G M; then there exists a unique nonnegative function Fa G 
HiociBa) n L^^J^{Ba; Ivfdri) solution of (2.6 ) and satisfying 

lim F„(j?) = oo. (2.7) 

Furthermore a Fa is decreasing. 

Proof Set K{r]) = el''!'/'*. Then (2.6 ) becomes 

- K-^div{KVf) - 7/ + I'?!'' \fr' / = 0. (2.8) 
Step 1- Boundary behaviour. First we claim that 

Actually, \i Q < h < \r}\ < a, u satisfies 

-K-^div{KVFa) - ^Fa + CFP < 
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with C = min{a'^, 6''}. We perform a standard variant of the two-sides estimate method used 
in [17] : we set T := Bp \ with b < p < a, a = {p — b)/2 and denote by z the solution of 



z" -CzP = in (-a, a) 
z{—a) = z{a) = DO. 



(2.10) 



Then z is an even function and is computed by the formula 

r , '^^ ^ = J^{a-t) Vte[0,a). (2.11) 

Notice also that lima^o z{t) = oo, uniformly on (—a, a) and 

,5,(,-„)V<.-.).«,^(J|±i>,)"""". (2.2, 

We set Z{t]) = z{\r]\ — {p+b)/2) and we look for a super-solution in F under the form w = MZ{r}) 
{M > 1). Then 

-K-'^div{KVw) - 7W + CwP = M (^{Mp-^ - 1)CzP - (^^-^ + ^) ~ ' 

Since 



^'w = - < c*z^p+^y\t), with c* = J^, 



we derive 

- K-^div{KVw) -'yw + CwP>M (^{MP'^ - 1)CzP - (^^^^ + |) C*z^p+^^/^ - 'yz^ (2.13) 

on {77 : (p — b)/2 \ri\ < p}; and the same inequality holds true on {77 : p < |ry| < {p — b)/2}, up 
to interverting a and 6. For any M > 1, we can choose 6 > such that for any b < p < a, the 
right-hand side of (2.13 ) is positive and maximum principle applies in Bp \ Bh. Thus MZ > Fa in 
r. Furthermore, the previous comparison still holds if we take p = a, which implies a = {a — b)/2. 
Therefore, using the explicit value of C 

Because M > 1 and < 6 < a are arbitrary, we derive 

limsup(aH.I)^/^-^)i^«(.)<(^4^ . (2.15) 

For the estimate from below we notice that u satisfies 

-K-^div{K\/Fa) - ^Fa + CFP>0 

in [i] : b < \ri\ < a}, with C = maxja'^, 6^^}. Taking now a = a — 6, we denote by z the positive 
solution of 

z" + jS-CzP = in (0, a) 

5(0) = (2.16) 
z{a) = 00. 
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Then z is computed by the formula 



POO 

Jz{t) 



ds 



^z'^iO) - + 2(7sP+i /{p + 1) 



a-t Vie[0,a), 



(2.17) 



and formula (2.12 ) is valid provided C be replaced by C. We fix A e dBa with coordinates 
(a, 0, 0), and look for a subsolution under the form w{r]) = Mz{r]i — h) with < M < 1. Then 



-K-^div{KVw) - 7w; + Cw^ = M ((M^-^ - 1)F 



since w' > 0. Applying again the maximum principle, we derive w{ri) < Fa in Bar\{r] : b < rii < a}. 
But clearly the direction rji is arbitrary and can be replaced by any radial direction. Thus 



liminf(a- |77|)2/(P-^)F„(r/) > M 



2(P+1) 



i/(p-i) 



In turn, (2.18 ) implies 



liminf(a- |r/|)2/(f-i)F„(77) > 



max{ aP/^, 6^/3 1)2 

2{p+i) 



aP/3(p_ 1)2 



(2.18) 



(2.19) 



and (2.9 ) follows from (2.15 ) and (2.19 ). 

Step 2- Uniqueness. If F' is another nonncgative solution of (2.6 ) satisfying the same boundary 
blow-up conditions, then for any e > 0, F^ = (1 + e)F' is a super solution. Thus, for 6 > 0, 

pp F'P 



div{KVFa) , div{KVFl) , , .p 

' ^ m 



Fa +5 

By monotonicity 



f; + s 

Fg ^ 

Fg + 6 F^ + S 



F' 



Fa + 5 F^ + S 
{{Fa + 5f-{F',+5f)+Kdn. 



K\{{Fa + 5f-{Fl + 5f)+dn 



FP 



and 



< 



Fa 



Fa+5 F^ 

f: 



F'P \ 

i{Fg+5f-{Fi + 5r)+>0, 



Fa + 6 F' + 6 



{{Fa + 5Y - (f; + s)')+ < {{Fa + sr - (f; + 5)% 



By Lebesgue's theorem, since (2.9 ) implies that {{Fa + 6)^ — {F^ + 5)^)+ has compact support in 

Ba, 

^" ^' ^ {{Fa + Sf-{F^ + 6y)+Kdv = 0. 



lim 

5^0 



Fa + 6 F' + 6 



Using Green formula, we obtain 

div{KVFa) div{KVF^) 



Fa + S 



F' + S 



{{Fa + 5f-{F',+5f)+Kd7i 



II 



VFa 



Fa>F' 



Fa + 5 

F' + 5 



VF' 



+ 



VF' 



F[ + S_ 

Fa + S 



Kdr] > 0. 
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Letting 6 ^ 0, we derive, by Fatou's theorem, 



7 — ^ ~ 7 

/ / {Fr' - K"-') {Fl - F^')Kdv < 0. 



Thus Fa < F^. Since e is arbitrary. Fa < F' . The reverse inequahty is the same. The monotonicity 
of a F(j is proved in a similar way, by the previous form of maximum principle. 

Step 3- Existence with finite boundary value. We shall first prove the existence of a positive 
solution Wk of (2.6 ) with boundary value equal to fc > for small value of a, and we shall let 
fc — > 00 in order to obtain one solution satisfying (2.7 ). We denote by Ja the functional defined 
over H^{Ba) n LP+\Ba; Ir/pd??) by 

JaH = (\yw\' - -fw^ + -l-\rjf\wf+'^ K{r])drj. 

Let A: > and k G C^{Ba) with < K{ri) < k, supp{k) (Z Ba\ Ba/2, = k on Ba \ -620/3- If 
V e H^{Ba) n LP+'^{Ba; \vfdr]) and w:=v + k, then 

Jaiw) = Ja{v + k)> Ja{v) + Ja{n) + I {Vv.Vk - 7m - \'nf\v\PK) K 

JBa 

Since 7 < Ao, it follows from Cauchy-Schwarz and Holder- Young inequalities that 

Ja{w) > {l-e^)Ja{v)-^Ja{K) 

for < e < 1. Because lima^a \a = 00, there exists ao G (0, 00] such that, for any < a < ao, 
Ja{v) is bounded from below on HQ{Ba)r\LP'^^{Ba; \r]\^dri). Thus there exists a minimizer Wk such 
that Wk = v + K with v in the above space; Wk is a solution of (2.6 ) and Wk\dB^ = k. Furthermore 
Wk is positive. Notice that if 7 < 0, oq = 00, in which case there exists a solution Wk for any fc > 
and any a > 0. The uniqueness of Wfe > 0, is a consequence of the monotonicity of the mapping 
k^ Wk that we prove by a similar argument as in Step 2: if A; < fc', there holds 

// - <r'){wl - wD) \nfKdn < o, 

which implies Wk <Wk- Uniqueness and radiality follows immediately, thus Wk solves the differen- 
tial equation 



—w — 



f N — lr\ 
( — h - 1 - 7W + r^wP = on (0, a) 



w{a) = fc and e Hl^,{Ba) n LP+l{Ba; Ir/l'^dT?). 



(2.20) 



Next we shall assume 7 > 0, cquivalcntly /? > —2. If Wk is a positive solution of (2.20 ) and A > 1 
(resp. A < 1) \wk is a super-solution (resp. a sub-solution) larger (resp. smaller) than Wk- Note 
that /3 > — 2 implies Wk{Q) > while /3 > — 1 implies also Wj.(0) = 0. Thus, by [13], there exists 
a solution w\k with boundary data Afc, and this solution is positive because Wk < wxk < Xwk 
(resp. Xwk < wxk < Wk)- Consequently, the set A of the positive d such that there exists a 
positive solution of (2.20 ) on (0, a) for any a < a is not empty and independent of fc. Furthermore, 
if for some a > and some fco > 0, there exists some positive Wkg solution of (2.20 ) on 0, a), 
then for any < a < d and any fc > 0, there exists a positive solution Wk of (2.20 ). Since 
r I— > max{fc, (7_|_a~'^)^/^^'~^)} is a super-solution, there holds 

Wk{r) < max{fc, {j+a-'^)^/^P-^^ Vr G [0,a]. (2.21) 
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Let us assume that a* = sup A < oo. Because of (2.21 ) and local regularity of solutions of elliptic 
equations, for any e, e' > 0, w^(a) is bounded uniformly with respect if e < a < a* — e'. But since 
(2.20 ) implies 



(a) = e^'-h^'^^w'^ie) + / {r^wl - ^wuV""-^ e^" 'Hr, 



w'i^{a) is actually imiformly boimdcd on [e, a*). It follows from the local existence and imiqucncss 
theorem that there exists 5 > 0, independent of a < a* such that there exists a unique solution z 
defined on [a, a + 5] to 

—— + -] z' -jz + rf^zP = on(0,a) ^2.22) 

z{a) = fc, z'(a) = w[,{o.), 

and S and k > can be chosen such that z > in [a, a + J]. This leads to the existence of a positive 
solution to (2.20 ) on [0,a + 5]. U a* — a < 6, which contradicts the maximality of a* . Therefore 
o* = 00. 

Step 4- End of the proof. We have already seen that fc 1— > Wfc is increasing. By Step 1, we know 
that, for any a > 0, and some b < a, there holds 

Wk{\v\) < C{a - |j?|)-V(p-i)on \ B,. (2.23) 

In particular 

Wk{b)<C* =C*{a,b,p,N) 

Next 

Wk{r) < max{C*, {-f+b-^Y^^P-'^^ Vr e [0,6]. (2.24) 

Combining (2.23 ) and (2.24 ) implies that w/. is locally uniformly bounded on [0, a). Since k ^ 
is increasing, the existence of Fa := i«oo = hm^^oo w'fc follows. The fact that a'l-^ decreases is 
a consequence of the fact that Fa' is finite on dBa for any a < a' . □ 

Remark. In the sequel we set F^o = lima^oo Fa- Then F^o is a nondecreasing, nonnegative solution 
of (2.6 ). Using asymptotic analysis, is is easy to prove that there holds: 

(i) ii/Sj^O 

Poo{v)=[j:^j |r?|-'5/(^'-i)(l + o(l)) asH^oo; (2.25) 

(ii) if /3 = 0, 

Fo.iv) ^ [j^j ■ (2.26) 

Furthermore, if /3 > —2, it follows by the strict maximum principle that Fa{0) = mm{Fa{r]) : 
\r]\ < a} > 0. This observation plays a fundamental role for obtaining estimate from above. 

Proposition 2.2 Assume p> 1 and (3 > —2. Then any solution u of (2.1 ) in Qt which verifies 

Vmiu{x,t)=Q Vx^O, (2.27) 

satisfies 

mx,t)| <min{c*|a;r('+'')/(^-'^i-('+'^)/2(p-i)j,^(^/^)| v(a;, i) e Qt \ {0}, (2.28) 
where c* = c* {N, p, j3) . 
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Proof. Let e > and a > and "Pa.c = {{x,t) : t > e, \x\/y/t — e < a}. By the previous remark 
minFa > 0, thus the function W{x,t) ^ (t - e)-(2+/3)/2(p-i)Fa(|a;|/Vt^), which is a solution of 
(2.1 ) in '?a,e tends to infinity on the boundary on 'Pa,e', since u is finite in Qt n ^a.e, W dominates 
u in this domain. Letting successively e and a ^ oo yields to m < F^o- The estimate from 
below is similar. Next we consider x G \ {0}, then v = \u\ satisfies (by Kato's inequality) 

dtv -Av + C{x)vP < in Bi^i/2{x) x (0, T), 

where C{x) = max{(|a;|/2)^; (3|a;|/2)''}. It is easy to construct a function under the form iu{y) = 
A (Ixp - 4:\x - yp)"^^^^""^'' which satisfies 

-Aw + C{x)wP = in B\^\/2{x) 

^^'^\x-y\^\x\/2W = 00, 

with A = A(a;) c* \x\^'^~'^^^^p^^\ c* = c*{N,p,(3) > 0. Using (2.27 ), it follows from Lebesgue's 
theorem that u{y,t) < w{y) in -B|^|/2(a;) x [0,r), thus u{x,t) < w{x) = c*|a;|-(2+/3)/(J'-i). Estimate 
from below is similar. □ 

The construction of the first part of the proof of Proposition 2.2 (estimate in To.e) shows that, 
without condition (2.27 ), equation (2.1 ) admits a maximal solution um- 

Proposition 2.3 Assume p > 1 and p > —2. Then any solution u to (2.1 ) satisfies 

|u(x,i)| < WM(a;,t) :=t-(2+/3)/2(f-i)Foo(a;M) V(a;, i) e Qt \ {0}. (2.29) 

As a variant of (2.28 ), we have the following Keller-Osserman type parabolic estimate which 
extends the classical one due to Brezis and Friedman in the case (3 = (see [3]). 

Proposition 2.4 Under the assumptions of Proposition 2.2 there holds 



c 

|a;|^+t)(i+/3)/2(p-i) 



< ^^-tt^t^^vt:;?::^;^ v(x,t) e Qt \ {o}, (2.30) 



with c = c{N, p, (3) . 



Proof. Assume \x\^ < t, then 



1 > 2-(2+/3)/2(p-l)f-(2+/3)/2(p-l) 



k|'+t)(2+/?)/2(p-l) 

2-(2+/3)/2(p-l) (2.31) 

mm{F^{r]) : |r?| < 1} 



2 

Assume |a;| > t, then 



> 2-(2+/3)/(p-i)|a;|-(2+/3)/(p-i). (2.32) 



(|a;|2+i)(2+/3)/2(p-l) 

Combining (2.31 ) and (2.32 ) gives (2.30 ). □ 



11 



2.2 Isolated singularities and the very singular solution 

Theorem 2.5 Assume p > 1 and —2<(3< N{p — 1) — 2. Then any solution u to (2.1 ) which 
satisfies (2.27 ) is identically 0. 

Proof. If -(2 + (3)/(p - I) + N - I > -1, equivalently (3 < N{p - 1) - 2, the function x ^ 
|2,|-(2+/3)/(p-i) jg locally integrable in M^, thus u{.,t) ^ in Lj^^{R^) as t ^ 0. For e > there 
exists R = R{e) such that u{x, t) < e for any |a;| > R and t > 0. Thus 

u{x,t + T) <e + E[uxB^u{;r)]{x,t) Vt > 0, r > and a; e M^, (2.33) 

where E[^] denotes the heat potential of the measure <p (see (2.2 )). Letting successively t ^ 
and e ^ 0, yields to u < 0. In the same way u>0. In the case /? = N{p — 1) — 2 estimate (2.30 ) 

reads 

\u{x,t)\< 



(|.T|'+i)^/2' 

From this estimate, the proof of [3, Th 2, Steps 5, 6] applies and we recall briefly the steps 

(i) By choosing positive test functions 0„ which vanish in V„ = {{x,t) : \x\^ +t< n~^} and are 
constant on = {{x,t) : \xf + t> 2n~^}, we first prove that, for any p > 0, 

{\u{x,t)\ + \xf\u\P) dxdt < 00. (2.34) 

SpX(0,T) 

Thus, using the same test function, we derive that the identity 

(^-udtC - uAC + \xf \u\^'^ uCj dx dt = 0, (2.35) 



i i Qt 



holds for any C e Co'^(M^ x [0, T)). The uniqueness yields to u = 0. □ 

Proof of Theorem, A- case I. In the case — 2 < /? < N{p — 1) — 2, the result is a consequence of 
Theorem 2.5. Next we assume /? < —2. If / is a solution of (1.6 ), it satisfies 

f{r]) = o(|77|-(2+/3)/(p-i)) as |?7| ^ oo. 
If /? = — 2, the equation becomes 

-A/-lr;.V/ + H-'|/rV = 0, 

and /(ry) ^ at infinity. Since any positive constant is a supersolution, / < 0. Similarly f > 0. 
If /3 < -2, for e > the function jj ^ e|r?|-(2+/5)/(p-i) = ^'(r/) belongs to Wl;^{R^) since /3 < -2 
and satisfies 



p — 1 J \p — 1 

Therefore, either if N >2 ov N = 1 and < — (p+ 1), V' is a super-solution of (1.6 ) for any e > 0. 

The conclusion follows as above. 

Finally we treat the case N = 1 and — (p+l)</3<— 2 where there exists a particular solution of 

r + ^/'+^^/-r^i/rv=o onM+, 



12 



under the form /i(r) = A/^j^pV (2+/3)/(p i)_ Furthermore, if / > (which can be always assumed 
by the maximum principle), it is a subsolution of the linear equation 

Noticing that this equation has a solution which has the same behaviour at infinity than the 
explicit solution of (1.4 ), namely 

<^i(r) = cr-(2+/?)/(p-i)(i + o(l)), 
by standard methods (see e.g. [10, Prop Al]), the second solution (p2 behaves in the following way 

02(r) = cr(2+/3)/(p-i)-ie-'''/4(l + o(l)) as r ^ oo. 

Consequently, by the maximum principle, any solution / of (1.4 ) on M such that /(r) = o((^i(r)) 
at infinity, verifies 

|/(r)| < C|rp+'^)/(f-i)-^e-'''/* for |r| > 1. (2.36) 
Using the equation, we obtain that 

fir) = e^'/'f (s^imr'm - ^/(«)) ds, 

thus 

|/'(r)| < Cr(2+/3)/(p-i)-2e-'-'/4 for \r\ > 1. (2.37) 
Since / e Hl^^{R), we derive that for any n e N*, 

/" (■^" " ^^') - - /(-")/'(-")) • 

Because of (2.36 ) and (2.36 ), this last term tends to as n — > oo. Therefore 

JZ " ^^') ''^^'"'^ = « ^ / = 0' 
which end the proof. □ 

Remark. The method of proof used in the case = 1 and —p — l</3<— 2is actually valid in 
any dimension, for any /3 < — 2. But it relies strongly on the fact that / e Hf^^, while the other 
methods use only / e W/„;^(R^). 

Proposition 2.6 Assume (3 > max{A^(p— 1) — 2; — A''} . Then for any k > there exists a unique 
solution Uk of (2.1 ) with initial data kdo- Furthermore k ^ Uk is increasing anduoo '■= limfe_>oo Uk 
satisfies Uoo{x,t) = t~^'^~^^^^^^P~^^ foo{x/y/i), where foo is positive, radially symmetric and satisfies 

f -A/oo - ir?.V/oo - 7/oo + \vf = in ,^ 
\ lim|,|^^|r?|(2+/3)/(p-i)/<^(ry)=0. 

Proof. The existence of Uk and the monotonicity of A; i— > Ufe has already been seen. By the uniform 
continuity of the Uk in any compact subset of Qt \ {(0, 0)}, the function Uoo satisfies 

liniUoo(a;,t) = \/x^0. (2.39) 
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For ^ > and u is defined in Qoo, we set 

Tf\u]{x,t) :=^(2+,3)/2(p-i)y(^^^^^)^ (2.40) 

If u satisfies equation (2.1 ) in Qoo, Te[u] satisfies it too. Because of uniqueness 

Te[Uk] = U£(2 + l3)/2(p-l)-N/2k- (2-41) 

Using the continuity of m i-^- Ti[u] and the definition of Uoo, we can let /c ^ oo in (2.41 ) and derive 
(by taking £t = 1 and replacing t hy £), 

T£[tioo] =Uo^=^ u^{x,t) = r(2+/3)/2(p-i)M^(x/VI,l). (2.42) 

Setting faoiv) = Uooix/Vt, 1) with rj = x/\/t, it is straightforward that f^o satisfies (2.38 ) (using 
in particular 2.39 ). Furthermore is radial and positive as the are. □ 

Lemma 2.7 The function foo satisfies 

fM = clTyl'^-^e-l^l'/^ (1 + o{\r,\-^)) as \r,\ - oo, (2.43) 

for some c = cn,p,i3 > 0. Furthermore 

.f^(^) = -^c|r;|27+i-^e-H'/4(l + o(|r;|-2)) as |r,| ^ oo. (2.44) 



Proof. Set r = I??! and denote foo{v) = /oo(^)- Then /oo satisfies, 

+ + +7/00 -r'^l/oorVoo = on (0,00), 

and linir^oo r"^^ foo{r) = 0. We consider the auxiliary equation 



(2.45) 



r+( ^^ + ^')/' + 7/ = on (0,00). (2.46) 



By [10, Prop Al], (2.46 ) admits two linearly independent solutions defined on (0,oo), yi and 2/2 
such that 

yi(r) =r-27(i + o(l)) and y2(r) =r2T-^e-'-'/4(i + o(l)), (2.47) 

as r — > 00. Next we choose R> Q large enough so that the maximaum principle applies for equation 
(2.46 ) on [i?, 00) and the yj arc positive on the same interval. For (5 > 0, Y5 = 5yi+foo{R)y2/y2{R) 
is a supersolution for (2.45 ). Furthermore fao{r) = o{Ys) at infinity. Letting 5^0 yields to 

fooir) < ^^y2{r) Vr > i?. (2.48) 

Using (2.47 ) we derive 

o<foo{v)<c\vf^-''e-MV^ y\n\>i. 

Plugging this estimate into (2.45 ), wc derive (2.43 ) from standard perturbation theory for second 
order linear differential equation [2, p. 132-133]. Finally, (2.44 ) follows directly from (2.43 ) and 
(2.45 ). □ 

An alternative proof of the existence of /oo is linked to calculus of variations. In the case (3 = 0, 
this was performed by Escobedo and Kavian [8]. This construction is based upon the study of the 
following functional 

Av) = U (\Vv\'-jv' + ^\vf\vf+')K{v)dv, (2.49) 
defined over the functions in B\{^^) n V^^^p^{R^). 



14 



Proposition 2.8 Assume p > 1 and (3 > N{p — 1) — 2. Then there exists a positive function 
/oo G nLj'+^^^(R^) satisfying 

-Af^-^r].Vf^-'yf^ + \vff^=0 m M^. (2.50) 

We recall that the eigenvalues of —K~^div{KV .) are the Xk = {N + k)/2, with k £ N and 
the eigenspaces Hk are generated by D"(j) where (f){r]) = K~^{ri) = e~l''l and \a\ = k. It is 
straightforward to check that J is C'^. In order to apply Ekeland Lemma, we have just to prove 
that J is bounded from below in H}^{R^). As we shall see it later on, the proof is easy when 
13 < N{p - l)/2, and more difficult when /3 > N{p - l)/2. 

Lemma 2.9 For any v G i?^(M^), there holds 

-( (2N +\r]\^)v'^K{r])dr]< [ \Vv\^K{r])dr]. 
Proof. We borrow the proof to Escobedo and Kavian. Put w = v\fK. Then 



\fKVv = Vw — —T]. 



Hence 



Because 



there holds 



/ \Vv\'^K{r])dri= I ( \\7w\'^ - wVw.rj + ^w'^ \ri\'^] dr]. 

— / wVw.rjdr] = [ w'^dr/, 

/ \Vv\^K{r])dr]= [ ( \Vw\^ + + ^w^ IrjA dr]. 
Jmn Jrn \ 2 4/ 

This implies the formula. □ 

Lemma 2.10 Let p > 1 and (3 < N{p — l)/2. For any e > there exists C = C{e,p) > and 
R = R{e,p) > such that 

r r / r \ 2/p+l 

/ v^K{r])dr]<e \VvfK{ri)dr] + C{ / \v\P+^\jjfK{r])] 
Proof. For R> there holds 

/ \ 2/(P+l) / X (P-1)/(P+1) 

/ v^K{'n)dr]<{ \v\P+^\r]fK{'n)dr]] i \r]\-^^/^P-^^ K{r])d'n] 

J\ri\<R \Av\<R J \Ari\<R J 

Since /3 < N{p - l)/2 N > 2(3 /{p - 1), we obtain 

(p-i)/(p+i) 



By Lemma 2.9 



/ |r,|-2/3/(f-i)if(,y)d,,) =C{R,N,p). 

J\r,\<R ) 

v^K{v)dri \VvfKirj)dri. 

7]\>R ^ Jr^ 
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The estimate follows by taking e = 4i? ^. □ 

It follows from the previous Lemmas that J is bounded from below in the space ff^(R^) fl 
whenever N{p - l)/2 - 2 < /3 < N{p - l)/2. Next we consider the case /3 > and we 
shall restrict the study to radial functions. 

Lemma 2.11 Assume /3 > 0. The functional J is bounded from below on the set 

X = e H}^{R^)r\ LP+pj^{R^) ■.v>0,v radial and decreasing | . 
Proof. For < ^ < i?, we write J{v) = Js,r{v) + J'g ^(w) + J'sniv) where 

l^yf - 7t;2 + -1- \rjf \vf+'^ K{r])dv, 



and 



^J5<\ri\< 

J'Uiv) = \l - + ^ \r,f \vr') K{rj)dv. 



Using Lemma 2.10, we fix R large enough so that J^'^ is bounded from below in H}^{R^) fl 
Lj'^^^(IR^). By Holder's inequality J'g j^ is bounded from below, thus we are left with Js,r- We 
assume that v is positive, radial, nonincreasing and v{S) = c = min{w(a;) : |a;| < S}. Then 

|^|P+i = ^P+i = (v-c + c)P+^ >iv- c)f +1 + and < 2{v - of + 2c^ , 
Js,R >U (\V{v - c)f - 2^{v -cf + ^ \rif \v - 0^+^) K{ri)drt + L{c), 

where 

L{c) = ^f \vf K{r])drj - [ K{rj)dr]. 

Clearly L(c)> M for some M independent of c. Therefore we are reduced to study the functional 
J5,R defined by 

JsA'^) = 11 f |V«^|' - 2jw' + \rjf \wf+'] K{n)dn 

over Hq j^{Bs) n L'^^^^ ^(Bg). Here we can fix (5 > small enough so that the first eigenvalue 

of —K~^div{KV .) is larger than 27, thus Js,r{v) is bounded from below in the class of radially 
symmetric nonincreasing, nonnegative functions v, and so is J. □ 

Lemma 2.12 Let v be a radially symmetric function in Hj^{W^)nL^^^^j^{W^). Then there exists 
a radially symmetric decreasing function v € Hj^{M.^) n L|'^^^(M-'^) such that J{v) < J{v). 

Proof. We define the two curves 

Ci = {{s,x) e M+ X R+ : -2-172;' + (p+ l)"'s^a;f+i = O} = {x = (2-1(^+1)75"'^)'/^^"'^} , 
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and 

C2 = {{s,x) e R+ X R+ : -72; + sf^xP = O} = {,x = (7.5-'^^''^^"''}. 

For fixed s > the function a; i-^- —2~^^x'^ + {p+ l)~^s^x^^-^ vanishes at x = 0. It has the foUowing 
properties: 

(i) it is decreasing for < a; < (75"'')^^'^ ^\ 

(ii) it achieves a minimum at Xg = {'ys~^)^^^^ ^\ 

(iii) and it is increasing for x > ['ys~^Y^^^ with infinite Umit. Furthermore it vanishes at 
Xg = {2-\p+l)js-^f^'-'\ 

Let w be a radially symmetric positive function. By approximation of radial elements in 
Hq j^{M^) n L|'^^^(M^), we can assume that v is with nondegenerate isolated extrema. 
We can also assume that the graph of v has at most a countable of intersections with C2, 
ai < 02 < as... < flfe < that the set of points {ak} is discrete, that all the intersections 
are transverse and that, for every j >0, 

w(s) < (7s '^j on (a2j,a2j+i), 

where ao = 0, and 

t;(s) > (7s '^j on (a2j+i,a2j+2+i)- 

The modifications of the function v is performed by local modification on each interval (afe, afc+i): 
Step 1- The construction of v on (02^,02^+1) is as follows. Let ai < a2 < ... be the sequence of 
local extrema of v, with v(ci2i+i) local minimum and i'(Q;2i+2) local maximum. By extension, since 
v'{a2j+{) > —^Idp— l)7^/'^~^^a^'^^^~^^^'^~^\ w(o2j+i) is a local maximum of v on (a2j,a2j+i)- 
If max{(a2i+i) : i > 1} < ^^(a2j+i), then v = max{w, u(a2j+i)}. 
If max{t;(a2i+i) : i > 1} > 'f^(ct2j+i), we define the increasing sequence {028^+1} by 

w(a2io+i) = niax{w(a2i+i) : i > 1}, 

w(a2ii+i)max{u(a2i+i) : % > io}, 

and by induction, 

v{a2ia+i)max{v{a2i+i) : i > id-i}- 

Thus we can assume that the local maxima of v are less than v{a2j+i) on the last interval 
{a2i^+i,a2j+i)- Next we define the function v hy v = max{v,v{a2io+i} on (02^, Q;2jo+i), v = 
max{?;,^;(Q;2ii+i} on (Q!2io+i) "2ii+i)- By induction, v = max{w, t;(Q!2id_i+i} on (Q!2i<j_i+i, Q!2i<j+i)- 
Finally v = max{w, i)(a2j+i)} on the last interval (a2id+i, 02^+1)- The function v is Lipschitz 

continuous, nonincreasing and, because v{s) < v{s) < {'ys~^)^^^^ ^\ there holds 



Ja2j<\v\<a2j+i V P+ ^ / 

< / (\Vvf-W + ^\vf\vr']K{rj)dv. 



(2.51) 



Step 2- The construction of v on (a2j+i, a2j+2) follows the same principle. Let /3i < /32 < ••• < (3d 
be the sequence of local minima of v on this interval. Furthermore i'(a2j+i) is the minimum of v 



on {a2j+i,a2j+2) and v'{a2j+2) < -P/{p - 1)7^^^^ 



(/3+p-l)/(p-l) 
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On (a2j+i,/3i) wc set v = mm{v,v{a2j+i)}- On (/3i,/32), v = min{v,v{(3i)}. By induction v = 
mm{v, v{(3i)} on On the last interval (/3d, ^2^+2), v = iiiin{w, z;(/3d)}. Because i) < 

on this interval and x 1— > —2~^'jx^ + {p+ 1)~^s^xP~^^ is increasing above the curve C2, we obtain 
similarly 



/ ( 

< 



,«~,2 



2 



/ 

"'023+1 



<l'/l<a2j+2 



■ ^v^ + 



p + 1 



(2.52) 



By construction v is nonincreasing. Combining (2.51 ) and (2.52 ), we obtain J{y) < J{v). 



□ 



Proof of Proposition 2.8. It follows from the previous lemmas that J is bounded from below on 
X and the function (p = belongs to X. Furthermore 



J{tcl>) = 



{N - 2j)t' 



jK-\n)d7i 



\t 



p+l 



Since /3>-/V(p— 1) — 2 «^=> N — 2^ < 0, the infimum m of J over radially symmetric functions 
is negative but finite and achieved by a decreasing function. Let {u„} C X a. sequence such that 
J{vn) I m. Then {vn} remains bounded in H^{M.^) n L'^^^^^{M.^). Up to a subsequence we can 

assume that w„ converges weakly in i/]^ (M^) and in _L|'^^a^(M^) and strongly in L]^(M^) to some 
function v. Moreover this convergence holds a.e., and, since Vn € X the same holds with v. Going 
to the limit in the functional yields to 



J{v) < liminf J{vn) — m; 



thus u is a critical point. 

The following uniqueness result holds. 



□ 



Proposition 2.13 Assume p > 1 and (3 > N{p — 1) — 2. Then foo = f<x- Furthermore foo is the 
unique positive solution of (2.38 ). 

Proof. We first prove that f^o is the unique positive radial solution of (2.50 ) belonging to 
H]^{M.^) n We denote r = |r/| and fooii) = foo{r). Let / be another solution in 

the same class. Thus there exists {r„} converging to 00 such that /(r„) 0. For e > 0, set 
/e = /oo + e- For n > no, large enough, «;_|_(r„) = 0, thus, as in the proof of Proposition 2.1, 



Vf-t^h 

Je 



V/. 



V/ 



Kdr] + 7 



if^-f^)+Kdr, 



We let successively r„ — * 00 with Fatou's lemma, and e 
e//e < 1 and (/^ - < + e L],{R^). We get 



+ // \vf ip-' - fr'){f' - m+Kdv < 0. 

with Lebesgue's theorem, since 



// 



/oo 



V/00 - y V/ 



+ \vf{f 



p—1 fP—^ 
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which impHes / < /oo- In the same way /oo < /• By Lemma 2.7, /oo e -ff^(M^) n ^[^^^^(IR^). 

Thus /oo = /oo- □ 

Wc end this section with a classification result 

Theorem 2.14 Assume p > 1 and /? > N{p — 1) — 2 and let u be a positive solution of (2.1 ) 
which satisfies (2.27 ). Then, 

(i) either there exists fc > such that u = Uk, 

(i) oru = Uoo- 

Proof. Because of (2.27 ) , the initial trace tr{u) of u is is a outer regular Borel measure concentrated 
at (see [12]). Then either the initial trace is a Radon measure, say kSo, and we get (i), or 



lim / u{x, t)dx = oo, (2.53) 

for every e > 0. This implies u > u^c. as in [11]. Notice that, in this article, this estimate is 
performed in the case /3 = 0, but the proof in the general case is the same. In order to prove that 
u < Uoo) we consider, for e > 0, the minimal solution v := of 

dtv -Av+ \xf\v\P-^v = in Qt 
tr{v) =Ub^, 

where ug^ is the outer regular Borel measure such that vg^ (E) = for any Borel set E c such 

that E O Be = 0, and lyg (E) — oo otherwhilc. This solution is constructed as the limit, when 
m ^ oo of the solution v^^m of (2.1 ) verifying Vc,m{-,0) = fnXg . Clearly u < v^. Furthermore, 
for any > 0, 

Te[Ve,m] = ^'e/VZ,m£(=+'3)/2(p-i) =^ = ^e/VI =^ ^^I^o] = ^^0, (2-55) 

where vq = limg^o'f^e- This, and the fact that \imt^oVo(x,t) = for every x G \ {0}, imply 
that vo{x,t) = t~^'^~^^^/^^P~^^ foo{x/Vi) = u^{x,t). At the end, since u < => u < vq, it follows 
u < Uoo- n 



3 Existence of very singular solutions 

Our study of the singularity set of the solution Uoo in the case of strongly degenerate potential 
(1.9) is based on some variant of the local energy estimate (abr. L.E.E.) method. First the L.E.E. 
method for the study of singular solutions of quasilinear parabolic equations was used in [15]. 
Adaption of this method to the study of conditions of removability of the point singularities of 
solutions of the quasilinear parabolic equations of diffusion-strong absorption type was given in [9] . 
In [14] there was elaborated a variant of the L.E.E. method, which allowed to find sharp conditions 
on the time dependent absorption potential, guaranteing existence of very singular solutions of the 
Cauchy problem to diffusion-strong absorption type equation with point singularity set. Here wc 
provide a new application of the L.E.E. method in describing the transformation of V.S.S solution 
into the R.B. solution in terms of the flatness of the absorption potential in the space variables. 
We consider the sequence of the Cauchy problems 

ut- l^u + h{\x\)\u\P-'^u = Q in x (0, T), p > 1, (3.1) 



u{x,Q) = UQ^k{x) = Mk exp(-2 Vo-^^fc)4(a;), 



(3.2) 
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where Sk is a regularized Dirac measure: 5k & C(M^), Sk S weakly in the sense of measures as 
fc — > oo, 

supp^fe C {x : \x\ < exp(-/iofc)} Vfc e N, (3.3) 
where the constant /xq > will be defined later on, and 

Mfe = expexpfc Vfc e N. (3.4) 
Without loss of generality we suppose that 

ll4|li,(Riv) < expifioNk). (3.5) 
We write the potential h in the equation (3.1) under the form, 

h{s) = exp(-a;(s)s-^) Vs > 0, (3.6) 
where u>{s) > is arbitrary nondecreasing function on [0, oo). 

Theorem 3.1 Let the function u){s) defined in (3.6) satisfy additionally the following Dini-like 
condition 

/ oj{s)s~^ds < ^2 < oo, di = const > 0, (3.7) 
Jo 

and the following technical condition 
sa;'(s) 



w(s) 



< 2 - ao Vs e (0, So), So > 0, < ao = const < 2 (3.8) 



Then the following a priori estimate of solutions Uk of the problem (3.1), (3.2), (3.5), holds uni- 
formly with respect to k Gf^, 



/ \ukix,t)\'- 



dx < Citexp 



C2($-M^ 



(3.9) 



where the constants Ci > 0, C2 > 0, C3 > do not depend on k. Here $ ^(s) is the inverse 
function to 

s^$(s):= / ^dT. 

r 



>(s) := /' 
Jo 

Let us define the following families of domains 

B{s):={x:\x\<s}, f2(s) := \ B(s), 



Qllis) :=n{s) X {ti,t2), Vs>0, \/0<ti<t2<T. 

Let u{x,t) = Uk{x,t) be a solution of the problem (3.1), (3.2) under consideration. We introduce 
the energy functions 

I{s,t):= f f {iV^ul"^ + h{\x\)\u\P+^)dxdt, (3.10) 
Jo Jn{s) 

J(s,t)= \u{x,t)\'^dx, E{s,t) = \uix,t)\'^dx. (3.11) 
Jn{s) Jb{s) 
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Lemma 3.2 The energy functions J{s,t), I{s,t) defined by (3.10), (3.11) corresponding to an 
arbitrary solution u = Uk oj problem (3.1), (3.2) satisfy the following a priori estimate 



_ (iV-l)(p-l) 2 

J{s,t)+I{s,t)<ctg{s):=ct[ / r p+3 h{r)p+^dr 



uniformly with respect to k (z N. 



P + 3 

p-i 



Vs > exp(-/iofc). (3.12) 



By c, Cj we denote different positive constants, which depend on known parameters N,p,aQ,d2 
only, but their value may change from lines to lines. 

Proof. Multiplying equation (3.1) by u and integrating in (5*J(s), we obtain the following starting 
relation after standard computations. 



\u{x, t2)\'^dx ■ 



n{s) 



{\V:,u\^ + h{\x\)\u\P+^)dxdt = 

= 2"^/ \u{x,ti)fdx+ f f u^dadt:=Ro + Ri. (3.13) 

Jn(s) Jti J\x\=s on 



Let us estimate Ri from above. Using Holder's and Young's inequalities we have 



/ u{x,t)p^da <cs^^^^^^ i [ \V,u\^da] if \uf+^ 

J\x\=s on \J\x\=s I V"'l^l='^ 



P+i 

'da < 



< cs%P+i) '\{s)-^ I [ (IV^wp + h{s)\u\P+^) da] 

\J\x\=s J 



p+3 
2{p+l) 



Integrating in t, we get 

du 
u— da dt 
J\x\^s on 



It is easy to see that 



(W-l)(p-l) ]^ p-1 

< CS 2(p+i) h{s) p-iT^I^+TT 



J\x\=s 



{\V^uf + h{s)\uf+'^)dadt 



p + 3 
2(p+l) 



(3.14) 



-^I{s,t)= r [ {\Vxu\^ + his)\u\P+') ds, -^J{s,t)>0. 

as Jo Ji^i=s as 

Therefore because of the property (3.3) satisfied by wo.fc, and estimate (3.14), we derive the following 
inequaUty from relation (3.13) with t2 = t, ti = 0, s > exp(— /xofc). 



p-l 1 (JV-l)(p-l) , ^ 

J{s,t)+I{s,t) < ct^(^h{s)~^s ( - — {I{s,t) + J{s,t)) 



p + 3 
2(p+l) 



(3.15) 



Solving this ordinary differential inequality (abr. O.D.I.) with respect to the function I{s,t) + 
J{s,t), we deduce that estimate (3.12) holds for arbitrary s > exp(— /UqA;). □ 

Next, we define Sfe > by the relation 



g{sk) = Ml" = exp(£o exp k), 



(3.16) 
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where < eq < 1 will be defined later on. Now we have to guarantee that 

Sfe > exp(-/iofc) := Sfe > ko{eo,ao,i'o,p). (3-17) 

Using [1, Lemma Al], it follows from the definitions (3.6)of function h{.) and (3.12) of function 
g{.), that the next estimate holds, 

^)"^9.is)<,is)<(^)'^M, (3.18) 

where gi{s) = s^~^~ p-i'a;(s)p^ exp ^ ^^^^^ '^o is constant from condition (3.8). The 

following simpler estimate follows from (3.18): 

for any s G (0,so), where sq = so(i^o) — > as j/q — > 0. As a consequence of definition (3.16) of Sk, 
and using (3.19), we get, 

— 2 7 7^ ^ ^0 exp k. (3.20) 

Sfe (P - 1) 

Integrating (3.8), we deduce that iv satisfies 

w(s) > s^-°"> Vs > 0. (3.21) 
Combining (3.21) and (3.20) we derive: 



ffo(p-l)/ V "0 

Next wc define from (3.2) and set /zq = 2ao ^. It follows from (3.22) that (3.17) is satisfied for 
all k > ko = ko{eo,ao,i'o,p). As result we derive that estimate (3.12) obtained in Lemma 3.2 is 
valid for s = Sk, i.e. 

J{sk,t) + I{sk,t) < ctg{sk) yk > k„ = kn{eo,ao,iyo,p)- (3-23) 

In order to find estimates characterizing the behaviour of the energy function E{sk,t) with 
respect to the variable t > 0, we introduce the nonnegative cut-off function (pk G C^(R) defined by 

V?fc(s) = 1 if s < Sk, ^pk{s) = Q if s>2sk, (Pk{s) < cs'^^ . (3.24) 

Multiplying (3.1) by Ufe(^^(|a;|) and integrating with respect to x, we get 

2"'^ / u\x,t)^l{\x\)dx + I \y,iu^k)\^dx+ I h{\x\)^l\u\P+Hx 

"I JRW J-g^N Js^N 

< [ u^{x,t)\V^ipk{\x\)fdx:=Ri. (3.25) 

By (3.24) and (3.23), we obtain 

Ki < cis^2 I \u{x,t)fdx < cisl'^J{sk,t) < C2SlHg{sk). (3.26) 

J Sk<\x\<2ak 
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Using (3.25), (3.26) and Poincare's inequality we derive the following differential inequality, 



di 



/ u^{x,t)(pldx] +dos^'^ u^{x,t)(pldx <-cs^'^tg{sk), do > 0. (3.27) 



We set 

^fe(i) := / \uk{x,t)f^pi{\x\)dx, 

and obtain the following O.D.I, from (3.27), 

V-feW + rfoSfe Vfc(*) < csftg{sk). (3.28) 
We rewrite (3.28) under the form 

i^kit) + Y^k^'^'^^t) + 2"' (dosfMt) - 2csftg{sk)) < 0. (3.29) 
Using the relations (3.2), (3.5) satisfied by Uk,o, we see that Vfe verifies, 

V'fc(O) < / \uk,o{x)fdx < Mk. (3.30) 

At last, we define the tk by 

tk = Msk) (3.31) 

where ui is the function in (3.6) and 7 > is a parameter which will be made precise in the next 
lemma. 

Lemma 3.3 There exists a constant 7 > 0, which does not depend on k, such that any solution 
^fc of problem (3.29), (3.30) satisfies the following a priori estimate 

i^kitk) < '^do^ctkg{sk) Vfc > k{so,i^o), (3.32) 
for some tk < tk, where tk is defined by (3.31). 

Proof. Let us assume that (3.32) is not true, and for any 7 > there exist k> ko such that 

Mt) > 2d^^ctg{sk) Vt:0<t< -fuj{sk) = tk- (3.33) 
This relation combined with (3.29) implies the following inequality, 

^kit) + Y'fMt)<0 yt:0<t<jw{sk). 
Solving this O.D.I, and using (3.30), wc get 



Mt) < V'fe(0)exp (-^) < Mke^p ("^) " ^'^^^'^^^ ^^'^^^ 
We derive easily the next estimate from (3.34) and (3.33) 

Mk exp l^ -^og(^fe) ^ > 2d^'cg{skhujisk). (3.35) 
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Using (3.16) and (3.4), we deduce from this last inequality, 

(1 - so) exp k > ^£2^^ + ln(2do ic7) - ln(a;(sfc)-i). (3.36) 

Similarly to (3.20), it follows, from (3.19) and the definition (3.16) of Sk, that there holds 

u>{sk) 2(1 + I'o) 



4 iP-i) 

Using this estimate and (3.36), we derive 



> eo exp k. (3.37) 



(1 - so) exp k > ''°7/f exp k + ln(rfo ^2c7) - ln(a;(sfc))-^ (3.38) 
4(1 + uo ) 

Noticing that (3.21) implies 

ln(w(sfc))-i < (2 - ao) ln(s^^), (3.39) 
and (3.22) can be writen under the form 

ln(s-)<^ln(f2i^V A, (3.40) 

^ ^ ao \2{l-i/o)J ao ^ ^ 

we deduce the following inequality from (3.39), (3.40) and (3.38), 



(2-ao)^ /eo(p-l) . ^^^^^^ 



(1 - eo) exp A; > ^"J/!^ , ^\^° expfc + ln(2do 'ci) - (2 - ao) — 
4(1 + I'o j ao 

2 

ao ""V2(l-i'o) 

If we define 7 by the equality 

^ X dojjp - l)£o (l-£o)(l + ;/o)8 _ 

(1 - £0) = on I ^ 7 = — J7 TT := 70, (3.42) 

8(1 + 1/0) do{p-l)eo 

then inequality (3.41) yields to 

(2-ao) , - /o^-i- ^ (2-ao), /£o(p-l)\ 

^ '-k > (1 - £0) exp fc + ln(2do C70) - ^ In — . 

ao ao V2(l-fo)/ 

It is clear that we can find fc = A:(£o, i^o) < 00 such that the last inequality becomes impossible for 
k > k, contradiction. Consequently, (3.33) does not hold for 7 = 70 and estimate (3.32) is true 
with 7 = 70. □ 

Proof of Theorem 3.1. Comparing definition (3.11) of E{s,t) and definition of V'fe, we easily see 
that 

E{su,t) < Mi) ^ E{sk,tu) < ^k(tk). (3.43) 
Therefore, using estimates (3.12), (3.32) and (3.43), we obtain 

/ \uk{x,tk)\'^dx = E{sk,tk) + J{sk,tk)<{do^c + c)tkg{sk). (3.44) 

Next we estimate the right-hand side of (3.44). Using (3.16), (3.31) and inequality (3.32), we get 

tkgisk) < 7ow(sfc)M^» < 7ow(so) exp(£oexpA:), (3.45) 
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where 70 is defined by (3.42)and sq > by (3.8). We obtain easily from (3.45) 

(cdo ^ + c)tkg{sk) < exp 
Let ki be tlie smallest integer such that 



^^^ln(7oa;(ao)(c + ccio-^)) ^^p^ 
exp K 



(3.46) 



ln(7oa;(so)(c + cdo^)) < eoexpfci, (3.47) 

equivalently 

fci = [ln(£oMn(7oa;(so)(c + crfo^)))] +1, 
where [a] denote integer part of a. Then it follows from (3.46) 

(cdfl ^ + c)tkg{sk) < exp(2£o exp fc) Vfc > fci. (3.48) 

If we fix £0 such that 

2£o < e-\ (3.49) 
then the next estimate follows from (3.44) and (3.45)-(3.49) 

/ \uk{x,tk)\^dx<Mk-i, (3.50) 

for all k > maxjfco, /c, /ci}, where /co is from (3.17 ), fc - from (3.32), and fci from (3.47). Estimate 
(3.50) is the final step of the first round of computations. For the second round, we begin by 
definiting Sk-i analogously to Sk'- 

g{sk-i) = Mll^ = exp(£o exp(fc - 1)). (3.51) 

Prom estimate (3.12) we obtain 

J{sk-ut) + I{sk-i,t) < ctg{sk-i), (3.52) 
since Sk-i > Sfe- Analogously to ipk, we define the function <fk-i and set 

ipk-i{t) ■■= / \uk{x,t)\'^\ipk-i{x)\'^dx. 

In the same way as (3.28), the following O.D.I, follows 

V'fc-i W + iV'fc-i(i) < cs^^M^k-i) > ife. (3.53) 

Using (3.50), we derive 

i'k-i(f^k) < Mk-i, tk<tk- (3.54) 

If we analyze the Cauchy problem (3.53), (3.54) similarly as problem (3.28), (3.30)) was analyzed 
in Lemma 3.3, we obtain the following a priori estimate for ipk-i{t), 

i^k-iitk +tk-i) < 2d^'^c{tk +tk-i)g{sk-i), (3.55) 
where tk-i < ife-i := 7ow(sfe-i), 70 is from (3.42). It is clear that 

E{sk-i,t)<^k-i{t) yt>tk, 
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consequently 

E{sk-i,tk + tk-i) <ipk-i{tk + tk-i) < 2dQ'^-c{tk + tk-i)g{sk-i)- (3.56) 
Prom (3.52 ), we deduce 

J{sk-i,tk+tk-i) + I{^k-i,tk + tk-i) < c{tk+tk-i)g{sk-i). (3.57) 
Summing estimates (3.56) and (3.57) we obtain 

\uk{x,tk + tk-i)\'^dx < {'cdQ^ + c){tk + tk-i)g{sk-i), (3.58) 

and we use this last estimate for performing a similar third round of computations. Iterating this 
process j times, we deduce 



Uk 



/ k-j 
K i=k 



/k-j 



dx < {cdg ^ + c) Xl^M 9{sk-j)- 



In particular, we can take j = k — I, where I G N satisfies 

l> lo := maxjfco, k, ki}. 



Then we obtain: 



i=k 



dx < {cd^ ^ + c) i^tij g{si). 



^i=k 



Next, we have to estimate from above the sum of the ti for which there holds 

I I 
^ti< ^7ow(si), 

i=k i=k 

where Sj is defined by g{si) = M?°. By the same way as in (3.37), we obtain 

2 ^ 2{1 + Uo)oj{s,) 2(l + i^V(so) , 

s| < — -, -T exp(-«) < — exp(-^) > Iq, 

(P - 1)60 (p - l)eo 

where lo is the integer appearing in (3.60 ), and from this inequality follows 

'2(l + z/o)w(so)^'^' 



Si < 



exp ( "2 ' Ciexp 



(P - 1)^0 

Therefore, using the monotonicity of the function w, we derive 



^^w(si) < w ( G\ exp 

i=k i=k 



<-j LO (Ciexp (^-^^^ ds 
y~'^uj{y)dy 



< 2 



Ci exp(-i^) 



ICi oxp(-§) 
-Ci exp(-i^) 

< 2 / y-^w{y)dy 
Jo 



:= 2$(C7iexp(-^i) 



(3.59) 

(3.60) 
(3.61) 

(3.62) 



(3.63) 



(3.64) 
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As a consequence of (3.62) and (3.64), we get 



^ < ^ < 270$ ( Ci exp( 

i=k i=oo ^ 



i-i 



) ■■=Ti. 



(3.65) 



The Dini condition (3.7) implies that — > as Z — > oo. Next, we deduce from (3.61) that 

/ \uk{x,Ti)\^dx<C2Tig{si), C2=cdo^+c Vfc > / > Zq- 
Using the fact that si : g{si) = Mf° and (3.66), we derive 

/ |ufc(a;,Ti)|^dx < C2T;exp(£oexp/). 

Because (3.65) imphes 

expZ = eC?(^$-i (^^' 



we get the following inequality by plugging last relationship into (3.67): 



T, 



I \uk{x,Ti)\'^dx <C2Tieyi-p 
At last, combining last estimate with (3.68 ), we obtain 
/ \uk{x,t)\'^dx < C2texp 



V/ > lo. 



>0, 



which ends the proof. 

Example 3.4 Assume u;{s) = < ao < 2. Then 



i>(s) = / S 

10 



^-"'°ds = 



2-ao 

Consequently, estimate (3.9 ) reads as follows. 



$~^(s) = (2 - a)^-''o s^- 



/ \ukix,t)\- 



dx < Cifexp 



C2 



2-ao 



> 0. 



(3.66) 



(3.67) 



(3.68) 



□ 



4 Razor blades 

In this section we consider equation (1.1 ) with potential /i(|a;|) of the form (3.6 ) with the limiting 
function w(|a;|) := |a;p£(|a;|), namely, we study the equation 



dtu -Au + e-^(l^l) \uf~'^ u = 0, in X (0, oo) 



(4.1) 



where £ e C(]R ) is positive, nonincreasing and limr^o^('') = oo. Our main result is the following 
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Theorem 4.1 Assume p > 1 and £ satisfies 

lim inf \x\'^i{x) > 0. (4.2) 

X— >0 

Then the solution Uk of the problem (1.1 ), (1.3 ), exists for any > and := limk^^ is a 
solution of (4-1 ) in Qoo \ {0} x M+ with the following properties, 

\imuao{x,t) = Va; ^ and lim Uoou(x,t) = oo Vt > 0. (4.3) 

Furthermore t t-^ Uoo{x,t) is increasing and \imt^ooUoo{x,t) = U{x) for every x ^ Q where 
U = linife^oo Uk and Uk solves 

-AC/fc + e-^(")L/^ = fc<5o tnV'{R^). (4.4) 

Proof. By assumption (4.2 ), property (1.2 ) is fulfilled. Thus for fc > there exists u := Uk 
solution of (4.1 ), (1.3 ). Moreover, for any k > there exists a solution Uk of (4.4 ) (see [18]); the 
mapping fc i— > C/fc is increasing and U = limfc_>oo Uk exists, because of Keller-Osserman estimate. 
U is the minimal solution of 

- Ay + e-^(^ = inR^\{0}, (4.5) 

veriiying 

/ V{x)dx = oo Ve>0. (4.6) 

If we denote by U the maximal solution of (4.5 ), it is classical that U = lime—^o where 

-AUe + e~^(^) UP = in \ J3e 
lim|^l_, Ue{x) =00. 

Since any Uk is bounded from above by U, the local equicontinuity of the Uk in (5t\{(0, 0)} implies 
that Moo satisfies limt_>o ^00(2;, t) =0 for all x ^0. 

Step 1: Formation of the razor blade. The Case 1: 1< p < 1 + 2/N . For e > 0, e-^d^l^ < e"^'^) 
for |x| < e. Therefore 

- Au + e-^(') u > 0, inS, x(0,oo). (4.8) 

and u > Ve in Be X {0,T) where Ve solves 

dtVe - Ave + e-^(<^) \vef~^ Ve = in x (0, 00) 

Ve = inaBeX(0,oo) (4.9) 
Ve{x, 0) = oo(5o in B^, 

where the initial condition is to be understood in the sense limfc^oo ^^o- We put 

We{x,t) = e2/(p-i)e-^(^)/(f-i)^;,(ex,e2t). 

Then We = w is independent of e and solves 

dtw — Aw + \wf~^ w = in i?! X (0, 00) 

w = in9Bix(0,oo) (4.10) 
w{x,0) = coSq in Si. 
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Therefore 

m(0, 1) > f.(0, 1) = e-2/(P-i)e^(^)/(f-i)w(0, (4.11) 
The longtime behaviour is given in [7] where it is proved 

lim e^'^w{0,T) = K(/)i(0). 

r— ^oo 

In this formula 0i is the first eigenfunction of —A in Wo'^(-Bi), Ai the corresponding eigenvalue 
and K> 0. Thus 

«(0, 1) > (5e-2/(P-i)e^(^)/(f-^)e^i^"Vi(0), (4.12) 
for some 5 > 0, if e is small enough. If we assume 



lim 



o\p-l 



2 , -1 ^(e) , _2 , 

Ine ^ + — ^ - Aie = oo, 



p-1 



(4.13) 



it implies 



«(0, 1) = oo u(0, t) = oo Vi>0. 



(4.14) 

Moreover, the unit ball Bi can be replaced by any ball Bji and Ai by Xn = R~^Xi. Therefore the 
suflacient condition for a Razor blade is that it exists some c > such that 



An equivalent condition is 



lini (^(e) — ce ^) = oo. 



liminfe^£(e) > 0. 



(4.15) 



(4.16) 



The gc;ii(}ial case. If p > 1 is arbitrary, we consider (3 > such that (3 > N{p— 1) — 2, and we write 

g-iix) ^ |^|/3g-^(a;)-/31n|a;|_ 

For R > small enough x i-^ £{x) :— £{x) + /3\n\x\ is positive, increasing and satisfies the same 
blow-up condition (4.2 ) as £. Clearly Uk is bounded from below on Bn x (0, oo) by the solution 
u := Uk of 

dtu -Au+ |a;|^e-^>) u = in Br x (0, oo) 

u = in 95;^ X (0, oo) (4.17) 
u{x,0) = kSo in Br. 

Therefore, for < e < i?, Woo is bounded from below on B^ x (0, oo) by the solution of 



dtVe - Ave + Ixfe-^'^"^ \vef \e = in X (0, oo) 

Ve = in dB^ X (0, oo) 
Ve{x,0) = oo6o in 



(4.18) 



If we set 



w,ix,t) = e(2+/3)/(f-i)e-^(^)/(P-i)^;,(ea;,e2i), 
then We = w is independent of e and 

dtw - Aw + Ixlf^ \wf~'^ w = inBix(0,oo) 
w = in dBi x (0, oo) 
w{x, 0) = oo(5o in Bi. 



(4.19) 
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By a straightforward adaptation of the result of [7] , there still holds 

lim e^i'"w(0,T) = k0i(O) 



for some k> 0. The remaining of the proof is the same as in case 1 < p < 1 + 2/N. 
Step 2: Asymptotic behaviour. A key observation is that, for any r > and any eo > 

/ Uoo{x, T)dx = oo. (4.20) 

We give the proof in the case 1 < p < 1 + 2/A/', the general case being similar. By step 1 

/ u{x,T)dx> [ t;,(a;,T)dx = e-2/(p-i)+^e^(^)/(f-i) /" w{y,e-^T)dy. (4.21) 

Jb, J Be JBi 



If we fix T and use [7], there exists eo such that w{y,e ^r) > 2 ^ne ^^'^ ^(piiv) ioi e < eo and 
yi G Bi. Therefore 

/ u(a;,r)da;>ce-2/(P-l)+^e^(^)/(^'-l)-^^^"^ (4.22) 
Jb, 

for some constant c > 0. If r is small enough, the right-hand side of (4.22 ) tends to infinity as 
e — > 0, so docs the left-hand side. This implies (4.20 ). For any fc > and any e > 0, there exists 
m = m(e) > such that 

min{u(a;, T),m}dx = k, 



L 



Be 

N 



thus, if we set (j)m = iiiin{u(a;, r), mlxs^ ) then u is bounded from below on M-™ x (r, oo) by the 
solution V = Vf^k of 



dtv -Av + e-^(^)|w|P-iw = in x (r, oo) 
v{x,t) = (pmix) in 



AT (4.23) 



When e — > 0, (f)m{.) kSo weakly in SDT(IR''^). By standard approximation property, v{e, k) — > Vo^k 
which is a solution of 

/ dtv -Av + e-^^'=^\v\P-^v = in x (r, oo) , . 

\ vi.,T) = k5o inR^. ^^-^^^ 

By uniqueness, vo^kix,t) = Uk{x,t — r). Letting fc ^ oo yields to 

u^{x,t + T)>Uoa{x,t) y{x,t)&QT. (4.25) 

This impHes that t Uoo{x,t) is increasing for every x € R''^. Because u{x,t) < U{x), it is 
straightforward that \imx^oou{x,t) — U{x) exists in \ {0}. 

Step 3: Identification of the limit. If C G C^(R^ \ {0}), there holds 

/ / f-u(a;,i)AC(a;)+e-^(^)uf(a;,t)C(a;))da;dt= / {u{x,T) - u{x,T +l))C,{x)dx. 
By Lebesgue's theorem 

/ f-C/(a;)AC(a;) + e-^(^)f/f (a;)C(a;)) dx = 0, (4.26) 
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and, from (4.20 ), 

U{x)dx = 00, (4.27) 



/ 

•> en 



for any eo > 0. Therefore [/ is a solution of tlie stationary equation (4.4 ) in M.^ \ {0} with a strong 
singularity at 0. For fc > and e > there exists k{e) > such that 



Jb. 



Uk(Adx = k. 



Let V := Vk,(i be the solution of 

f dtv- Av + e-^'^='">\v\P-^v^Q in Qt, 



piv (4.28) 



Since f/j ,.(.,0) < J7fe(g)(.), the maximum principle implies Vk.t < Uf^i^^y If we let e ^ 0, Vk,e 
converges to the solution Uk with initial data kSo- Furthermore k{e) ^ oo as e ^ 0. Therefore 

Uk{x,t)<U{x) 'i{x,t)eQT. (4.29) 

Letting successively fc — > oo and f — > oo implies 

U{x) < U{x) Vx e M^. (4.30) 

Since U is the minimal solution of (4.5 ) verifying (4.6 ), it follows that U = U. □ 
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